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ABSTRACT 

Global stability of an epidemic model for vector-borne disease was studied by Yang et al. [7]. A reinvestigation of 
the model with a saturated incidence rate and a treatment function proportionate to infectious population I is presented to 
understand the effect of the capacity for treatment. An equivalent system is obtained, which has two equilibriums; a 
disease-free equilibrium and an endemic equilibrium. The stability of these two equilibriums can be controlled by the basic 

reproduction number 9? 0 . The global stability of the disease-free equilibrium state is established by Lyapunov method and 
a geometric approach is used for the global stability of the endemic equilibrium state. The model has a globally 
asymptotically stable disease-free solution whenever the basic reproduction number is less than or equal unity and has 

a unique positive globally asymptotically stable endemic equilibrium whenever 5R 0 exceeds unity. Numerical examples 

are given for the model with different values of the parameters. Graphical presentations are also provided. The details are 
supplemented by numerical results given in annexure. 

KEYWORDS: Epidemic Model, Vector-Borne Disease, Saturated Incidence, Equilibrium Point, Stability, Reproduction 
Number, Treatment Function 

AMS Subject Classification: 92D25, 92D30. 

1. INTRODUCTION 

The main purpose of this paper is to study the dynamics of Vector-borne disease and understand the effect of the 
capacity for treatment. Vector-borne diseases are infectious diseases caused by viruses, bacteria, protozoa or rickettsia 
which are primarily transmitted by disease transmitting biological agents (anthropoids), called vectors, who carry the 
disease without getting it themselves. Globally, malaria is the most prevalent vector-borne disease whose vectors are the 
mosquitoes. The mosquitoes are vectors of a number of infectious diseases most prominent among which are dengue (the 
second most important vector-borne disease), yellow fever, St Louis Encephalitis, Japanese Encephalitic, and West Nile 
Fever, caused by the West Nile Virus. The literature dealing with the mathematical theory and dynamics of vector-borne 
diseases are quite extensive. Many mathematical models concerning the emergence and reemergence of the vector-host 
infectious disease have been proposed and analyzed in the literature [9, 10]. The mathematical models of epidemiology for 
the West Nile virus and for dengue have been investigated in [1, 3] and [10, 11]. 

Mathematical modeling became considerable important tool in the study of epidemiology because it helped us 
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to understand the observed epidemiological patterns, disease control and provide understanding of the underlying 
mechanisms which influence the spread of disease and may suggest control strategies. The model formulation and its 
simulation with parameter estimation allow us to test for sensitivity and comparison of conjunctures. Treatment plays an 
important role to control or decrease the spread of diseases. Kar et al. [16,17] and Ujjainkar et al. [4] proposed an epidemic 
model with treatment function. Wei et al. [6] investigated an epidemic model of a vector-borne disease with direct 
transmission and the vector-mediated transmission. The global stability of an epidemic model for vector-borne disease was 
studied by Yang et al. [7]. Here we have reanalyzed the model [7] with saturated incidence and a treatment function 
proportionate to infectious population I. The only mode of transmission in this model is through the vector. 

The paper is organized as follows: In Section 2, a vector-host epidemic model with vector transmissions is 
presented, where the dynamics of the hosts and vectors are described by SIR and SI model, respectively. Equilibrium points 
and reproduction number are obtained in Section 3. The analysis of stability of the equilibrium of the model is investigated 
in Section 4 and Section 5. Numerical analysis and conclusion are given in Section 6 and Section 7. 

2. THE MATHEMATICAL MODEL 

In this section, a vector-host epidemic model with vector transmissions is presented and investigated, where the 
dynamics of the human hosts and vectors are described by SIR and SI model, respectively. The total host population 

A^COis partitioned into three distinct epidemiological subclasses which are susceptible, infectious and recovered, with 
sizes denoted by 5(f) , /(f) and R(t) , respectively, and the total vector population /V, (f ) is divided into susceptible 
and infectious, with the sizes denoted by M(t ) and V (t) , respectively. 



The dynamics of this infectious disease in the host and vector populations are described by the following system 
of nonlinear differential equations: 

dS , 1SV 

— = b, ^ jU.S , 

dt 1 (1 + ay) 



dl _ \SV 



dt (1 + ay) 
dR_ 
dt 
dM 



-yl-jUj-rl, 



dt 



yl-^R + rl, 

■b 2 -A 2 MI -jU 2 M, 



( 2 . 1 ) 



dV 



With the initial Conditions S (0) = S 0 , I (0) = I Q , R(0 ) = R () , M (0) = M {] and V (0) = V Q . 

The parameters in the model stand for 

b ] = Recruitment rate constant in the host population 

b = Recruitment rate constant in the vector population 
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//, = Death rate constant in the host population 

//, = Death rate constant in the vector population 
y = Recovery rate constant in the host population 

\ = Transmission rate from infected vector population to susceptible host population 
/?2 = Transmission rate from infected host population to susceptible vector population 



<Tj =Constant parameter 



r = Treatment rate constant in the host population. 

It is assumed that all parameters are positive. The total dynamics of host population can be determined from the 
differential equation dN l / dt =b { — /il^ N ] which is derived by adding first three equations of system (2.1). The total 

dynamics of vector population can be determined from the differential equation dN 2 / dt = b 2 — jU n N 1 which is 
derived by adding last two equations of system (2.1). It is easily seen that both for the host population and for the vector 
population the corresponding total population sizes are asymptotically constant: lim N i = / fl x and 

lim N 2 = b 2 / fl 2 . This implies that in our model we can assume without loss of generality that N ] = b ] / jl x , 

t — »°o 

N 2 = b 2 t fl 2 for all t > 0 provided that 5(0) + 1 (0) + R( 0) = b { / fl x , M (0) + V(0) =b 1 / . 

Therefore, we attact system (2.1) by studying the subsystem given by 



dS 4 5V 

— = b \ M\S 

dt (1 + ay) 



dl_ 

dt 

dV 

dt 



\SV 
(1 + a x V) 



-yl-juJ-rl, 






C b 5 

— -V 

V^2 J 



I-M 2 V. 



From biological considerations, we study the system (2.2) in the closed set 

T = {(.S', /, V) eR 3 + : 0<S + / <bj //, , 0 < V <b 2 ! ju 2 ,S>0,I>0], 



(2.2) 



Where R denotes the nonnegative cone of R 3 including its lower dimensional faces. Obviously T is 
positively invariant set of (2.2).We denote by 3T and F° the boundary and the interior of T in R 3 respectively. 

3. MATHEMATICAL ANALYSIS OL THE MODEL 

Equilibrium points of model (2.2) can be obtained by equating right hand side to zero. The system (2.2) has two 
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equilibrium states: the disease-free equilibrium E {) = ( — - , 0, 0) G 3F and a unique endemic equilibrium 

A 

A = (s*,/\v*)e r° with 



5* 



(y+ /i { +r)(\ + ay*)f 

Ay* 



(3.1) 



A l A 2 b l b 2 -jU lJ u;(y+jU l + r) 

A 2 (Y+U+r)[A ] b 2 +u(n 2 + ay)Y 

A (a +^J ) 



(3.2) 



(3.3) 



The dynamics of the disease are described by the basic reproduction number 5K 0 .The threshold quantity 9? 0 is 

called the reproduction number, which is defined as the average number of secondary infections produced by an infected 
individual in a completely susceptible population. The basic reproduction number of model (2.2) is given by the expression 

9? 0 = . (3.4) 

AA(r+A + 0 

For < 1 , the only equilibrium is disease-free equilibrium E {) in r)T and for 9t () > 1, there is a unique endemic 

equilibrium E in T" . 

4. STABILITY OF THE DISEASE-FREE EQUILIBRIUM 

In this section we discuss the stability of the disease-free equilibrium EL . 



Theorem 4.1: The disease-free equilibrium E Q of (2.2) is locally asymptotically stable in T if^A Q <1; it is 
unstable if > 1 . 

Proof: To discuss the stability of the system (2.2) the variational matrix is 



J = 



A,v 



1 + ay 
Ay 



+ J u l 



1 + ay 

0 



0 

~(Y+ fi x + r) 
f b 



A X S 



(1 +ay) 

A,. S’ 



A 



vA 



(1 + ay) 

-v I -(a,/ + a) 



Z7 U \ 

At the equilibrium point = ( — , 0, 0), the variational matrix becomes 

A 
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"A 



J(E 0 ) = 



0 



0 

v 



0 

-(r + Mi + r) 
A 2 b 2 

Mi 



Mi 

*i b i 

Mi 

~ ft 2 



Its characteristic equation is 



(i+AP + (7+A + r + Mi)^ + Mi (7+Mi + f)~— -^—E] -0 

MiMi 

=> (^+A ) + (y+ A + r + Mi)^ + Mi(Y+ Mi + r )(l - ^o)] = 0 (4-1) 



By Descartes rule of sign all roots of equation (4.1) are negative if 9? 0 < 1 .Thus, if 9? 0 < 1 then the disease-free 
equilibrium E 0 is locally asymptotically stable; Otherwise, if 9? 0 > 1 then it is unstable. 

Theorem 4.2: If 1R 0 < 1, then the disease - free equilibrium E Q is globally asymptotically stable in T. 

Proof: To establish the global stability of the disease-free equilibrium, we construct the following Lyapunov 

function: 

L = M t MiJ + A ] b\V . 

Clearly, L> 0 along the solutions of the system (2.2) and is zero if and only if both I and V are zero. 
Calculating the time derivative of L along the solutions of system (2.2), we obtain 



L = M t Mi !+ \b x V 
^SV 



- M\Mi 



(1 + a,V) 



-{r+Mi+r)! 







fh, > 




1 

+ 

J3- 


K 


-2— y 


I-MiV 




\Mi J 





h 



<MiMA—y 

Mi 



MiMiir+Mi + r )- ^ 2 ^ 1 

Mi 

= -I [AA_ ( 7+ Mi + r) (1 - ' 9^o ) + ] 

< 0 , 



/ - \ypyi - \b { M,v 



1 b 

Where in the first inequality we have used the fact that < 1 and S < — in T. In addition, the last 

(1 + ay) Mi 

inequality follows from the assumption that S? 0 < 1 . Thus L(t ) is negative if 9t 0 < 1 .When 9? 0 < 1 , the derivative L = 
0 if and only if / =0, while in the case 9? 0 = 1 the derivative L = 0 if and only if / = 0 or V = 0. Consequently, the 
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largest compact invariant set in {(S', I, V) E T: L = 0} when 9t 0 < 1, is the singleton { E {) } . Hence, LaSalle’s invariance 
principle [8] implies that E {] is globally asymptotically stable in L. This completes the proof. 

STABILITY OF THE ENDEMIC EQUILIBRIUM 

In this section we discuss the stability of the endemic equilibrium E* . 

Theorem 5.1: If > 1, then the endemic equilibrium E of (2.2) is locally asymptotically stable in T 0 . 



Proof: At the equilibrium point E* , the variational matrix becomes 



/(£*) = 



Ay* 



v \+ay 



-+A 



\s* 



Ay* 



1 + a,V* 



-(Y+Mi+r) 



K 



r b ^ 

ZL-V* 

v Ei j 



(l +ay*Y 
A l S* 

(l +ay*) 2 
~{A i r* +// 2 ) 



Its characteristic equation is 

+ G\A~ + Cl 2 A + = 0 , 



(5.1) 



Where 



a. = 



Ay* 



1 + ay 



- + 2//| + y + r + /i 2 + A-,1 > 0, 



Ay* 



\ + a,V 



- + E, 



> 



(l + ay*) 

A,V 



(Y+E i+r) + 

* f u \ 

-v * 

'k.Ei j 



Ay* 



1 +ay 



-+2//j + y+ r 



U+v) 



w [ h r 



1 +ay* 



+ 2 //, + 7 + r 












1 +ay 



f b ^ 



+ A 



(/+//i+// 2 +r) 



Ml 






4^* 



l + ay ! 



- + 2//j + y+ r 



a 2 i* + 



Ay* 



1 +ay 



+Mi 



(Y+ff+Mi + r) 



>0 



[from (2.2)] 
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a 3 = 



xy' 



1 +ay 



+ M\ 



{y + Mi + r)(jU 2 +X_I*) 





( b \ 

2 V* 


(l + ay*) 


2 


v Mi J 



>(r+Mi + r) 



wr , W , .. ; r 

( 1 +ay*) (i+ay*) 1 2 



+ MiMi(Y+Mi + r)- 



MV 



(l + ay*) 



jU 2 



J 



={y+ J u 1 + r) 
>0 



Ay 2 v*r 






(l + a,r) (l + flfjV*) 



[/ro,77 (2.2)] 



It can be easily seen that a v a 2 ,a 2 > 0 . By Descartes rule of sign all roots of equation (5.1) are negative. Thus, 
the endemic equilibrium E is locally asymptotically stable in F° . 

Now, we analyze the global behavior of the endemic equilibrium E . Here, we use the geometrical approach of 

Li and Muldowney [14] to investigate the global stability of the endemic equilibrium E in the feasible region T. We have 
omitted the detailed introduction of this approach and we refer the interested readers to see [14]. We summarize this 
approach below. 

Let X f ( X ) E R" be a C* function for x in an open set I) c R" Consider the differential equation 

X=f(x). (5.2) 

Denote by x (t, X Q ) the solution of (5.2) such that x(0, X Q ) = X () . We have following assumptions: 

( ) D is simply connected; 

(H 2 ) There exists a compact absorbing set K c Z); 

( H 3 ) Equation (5.2) has unique equilibrium X in D. 

Let P : x \— > P(x) be a nonsingular ( " ) x ( " ) matrix-valued function which is C* in D and a vector norm I • I 
on R N , where N = ( n ). Let /i be the Lozinski“i measure with respect to the I • I. 



Define a quantity r/, as 
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q 2 = limsup sup -\ /u{B{x(s,x 0 )))ds, 

t—>oo x a —>K t 0 



(5.3) 



Where B = P f P 1 + PJ 1 2 P 1 , the matrix P f is obtained by replacing each entry p of P by its derivative in the 

direction of f ( p rj ) , and is the second additive compound matrix of the Jacobian matrix J of (5.2). The following 
result has been established in Li and Muldowney [14]. 

Theorem 5.2: Suppose that ( H ] ), ( H 2 ) and (H 3 ) hold, the unique endemic equilibrium E is globally stable 
in r° if q 2 < 0. 



Obviously F° is simply connected and E is unique endemic equilibrium for >1 in T 0 . To apply the 
result of the above theorem for global stability of endemic equilibrium E , we first state and prove the following result. 

Lemma 5.3: If > 1, then the system (2.2) is uniformly persistent in T 0 

System (2.2) is said to be uniformly persistent [2] if there exists a constant c > 0, independent of initial data in 
r°, such that any solution (5(f), 7(f), V(t)) of (2.2) satisfies 



liminf S >c , liminf I > c , liminf V >c 

t—>oo f— >00 ?— >00 



Provided (5(0), 7(0), V(0)) € T" . 

The uniform persistence of (2.2) can be proved by applying a uniform persistence result in [5, Theorem 4.3], and 
using a similar argument as in the proof of proposition 3.3 of [13]. The proof is omitted. 

The boundedness of F° and the above lemma imply that (2.2) has a compact absorbing set K c F° [2], Now we 
shall prove that the quantity q 2 < 0. 



/ \ 3/ 

Let v = (S,I, V ) and fix) denote the vector field of (2.2). The Jacobian matrix J = f— associated with a general 

OX 



solution x(t) of (2.2) is 



J 



+ /A 

) 



y 1 + CfjF 



W 

1 +ay 



0 



0 

-iy+Mi+r) 

fb ^ 

X, — — V 

J 



\ s 



(1 +ay)~ 

( 1 +ay ) 2 
-(Xj +M 2 ) 



and its second additive compound matrix J ' is, 
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r[2] _ 







f a 


1 + a 22 


a 23 


-a A 

U 13 








/ [ 2] = 




a 32 


j 1 


a \2 


9 








\ 


~ a 31 


a 2l 


a 22 a 33 y 






f Ay 


+ 2 ih + 7+ 


A 




As 




As 


\ 


[l+ojV 


j 




(i+ory) 2 




(i+ay) 2 





A 



( b 

-2L-V 

2 

0 



-( ^ V +//,+!,/ + //,) 
1+ay 1 - - 



i+<*y 



“(/+ A +f + /^/ + //,) 



Set the 



function P(x) = P(S, I ,V) = diag j 1, — , — 1 , then P,P 1 = cliag \ 0, , 1 . The matrix 

v ; \ v v\ f i v i v 



B = P f P + P.J - P can be written in block form as 



B 



B u B n 

\^2l ^22 J 



With B u =- 



( Ay } 

+ ay J 



^12 - 



\sv \sv 



(1 + ay) I {1 + ayj I 



&21 — 



\^2 J 

0 J 



V 



f i V ( xy . . A 

+ A + 'V + Bn 



B 



22 



I V 



1 + ay 
Ay 






\ + a,V 



— - — -(y+B ] + r + A 2 I+p 2 ) 



Consider the norm in R ' as I (u, v, w)l = maxflwl, Ivl + Iwl) where ( u , v, w) denotes the vector in R . Let g denote 



the LozinskiT measure with respect to this norm. Using the method of estimating ji in [15], we have 

fl{B) < sup(g p g 2 ), 

Wheie gj — // 10 (B n ) + |fi| 2 | , g 2 — |-® 2 i| 



(5.4) 
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1^12 1 ’ |^21 1 are matrix norms with respect to the Z 1 vector norm, and jU m denotes the Lozi nski“i measure with 
respect to the norm. 

From system (2.2) we can write 

/ \sv 



1 (l+ay)l 



ir+Mi + r), 



v_ 

V 



a 2 



( b \ 
\Ml 



V 



-Mi’ 



(5.5) 



(5.6) 



Since B n is a scalar, its Lozinski'i measure with respect to any vector norm in R will be equal to B n . Thus 



M\o(B\ i ) 



Av 



y 1 + tTjV 



+ 2//, + y+r 



in l \ SV 

Also | /}.. = i — . 

1Z I / - 



(1 + a { vy I 

On substituting the values of // |() ( B ] ( ) and |B p |, gj will become 



8i = 



< 



A,v 



y 1 + ay 



+ 2//, + y+r 



+ ■ 



A l SV 



(l + ay ) / 



4V 

i+ay 



+ 2// x + y+r 



+ 



A l SV 



(l +ay)i 



r Ay . W ' , 

+ 2 juyy+r h \-[y+ / j a +r) 

J 1 



y 1 + ay 



( 



\v 



- + /U, 
V 1 + ay J 



[From (5.5)] 



< 



I_ 

I 



Mi- 



In I V' b 
Now \B-,, = — — 

! y 



A 



— — V 

V Mi j 



(5.7) 



. To calculate /f |(| (/?,,) , we add the absolute value of the off- diagonal elements 
to the diagonal one in each column of B 21 , and then take the maximum of two sums, see [18], which leads to 



I v 

M W (B 22 ) = sup y- — - 



__yy_ 

v i+a,y 



+ Mi + A 2 I + Mi +-^ L —, — - — -(y+M l + r + \I+Mi) 1- 



1 + a,V I V 
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- su pjy~ — ~(/A +V + /0>y— -(r+A + r + V + /'2)i 



/ V , , . V 

J y ( /^2 / • 



On substituting the values of jU 10 (B 21 ) and S 21 , g 2 will becomes 



#2 



V 



— V 



jU 2 



+ — - — +A 2 I + // 2 ) 



V 1 V ( IT \ 

- — + Hi + — — y — (//, + Xj + // 2 ) 



[From (5.6)] 



-y~(/A + / V) 






( 5 . 8 ) 



Using equations (5.7) and (5.8) in equation (5.4), we have 



ju(B)<sup(g l ,g 2 )<j-jU l . 



(5.9) 



Along each solution x(t, X () ) to (2.2) such that X (} £ K , the absorbing set, we have 



-jju(B)dt < -log 
to t 1(0) ^ 



(5.10) 



Which further implies that 



- r l, m 

q 2 < lim sup sup - log — — - Hi 

t — >°° x a — >K t 1 (^vJJ 



= -A 

<0 

T-I* 

i.e. q 2 < 0. Therefore all the conditions of Theorem (5.2) are satisfied. Hence unique endemic equilibrium h is 
globally stable in F° . 

NUMERICAL SIMULATIONS 

In this section, we give numerical simulations supporting the theoretical findings. When we choose the values of 
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the parameters as: b { = 200 , b 2 = 300 , = 0 . 5 , /U 2 = 0 . 6 , Y — 0 . 4 , A^ = 0 . 01 , A^ = 0 . 02 , OC { = 0 . 1 , r = 0.4 for the 

model then £ (5 = 335 . 7049 ,/ = 24 . 72888 , V = 225 . 9217 ) exists and = 51. 28 > 1 . Our simulation shows 

endemic equilibrium E is asymptotically stable (see figure 1). To see the dependence of the steady state value / of the 
infective population on the parameter V, we have plotted figure 2 for different values of V, keeping all others parameter 
values same as for figure 1 and see that the infective population decreases as the parameter V increases. Further, we have 
also plotted figure 3, 4, 5 to see the dependence of steady state value of the susceptible and infective population on the 

parameter CC ^ , keeping all others parameter values same as for figure 1 and noted that the susceptible population increases 
and infective population decreases as increases. The details are supplemented by numerical results given in annexure. 



Stability graph for the point 
(S=335.7094, 1=24.72888, V=225.9217) 
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Figure 1 



Stability graph for I at different values of r 
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Figure 2 



Impact Factor (JCC): 2.0346 



NAAS Rating: 3.19 



Modeling and Analysis of a Vector-Host Epidemic Model with Saturated Incidence Rate Under Treatment 



13 



Stability graph for S at different values of a, 
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Figure 3 
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Figure 5 



CONCLUSIONS 

In this paper, we have studied a vector-host epidemic model with saturated incidence and a treatment function 
proportionate to infectious population I. The global stability of the disease-free equilibrium state is established by 
Lyapunov method and a geometric approach is used for the global stability of the endemic equilibrium state. The basic 
reproduction number is obtained and it completely determines the dynamics of the ODE model. The model has a globally 

asymptotically stable disease-free solution whenever the basic reproduction number 9^ 0 is less than or equal unity and has 
a unique positive globally asymptotically stable endemic equilibrium whenever 9? 0 exceeds unity. However, it is clear 

that when the disease is endemic, the steady state value I of the infective individuals decreases as the treatment function 
and d x increases, and I approaches zero as the treatment function and Ot x tends to infinity. Thus, it will be of great 
importance for public health management to maintain the treatment and saturation effects. In the absence of the treatment 
function and with ( X { = 0, the result is perfectly in agreement with Yang et al. [7]. 
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